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INTRODUCTION 



1.1 GENERAL. - - Marine navigation is the process 
of directing the movements of a vessel from one point to 
another. Historically, niloting, dead reckoning, and 
celestial methods have been used for this puroose. The 
earliest form of marine navigation was piloting, which 
came into being when man guided his vessel with resoect 
to familiar landmarks. Dead reckoning probably followed 
as he attempted to predict his future positions in venturing 
beyond sight of land. While it is generally accepted that 
celestial bodies were used to steer by from almost the 
beginning, celestial navigation, as it is known today, had 

to wait until information regarding the motion of celestial 

k e- • a [,J 

bodies was acquired. 

The advent of World War II brought about the development 
of various navigational aids and electronic oositioning 
systems which enabled the navigator to fix his position 
without regard to cloud cover or weather conditions. 'lince 
this time the trend in navigation eouinment has been toward 
obtaining greater positional accuracy, continous position- 
fixing capabilities, and lonp:er range operations. These 
improvements are being achieved by the use of electromagnetic 
and acoustic signals, ^oohisticated electronic and celestial 
methods, and inertial systems. Automatic computation has 



1 



' 

I 

i 

I 












become a necessity in the more sophisticated systems which 
require long and tedious calculations. 

The navigation requirements of the Navy are auite 
different from those of the ordinary user. an ideal system, 
other than a completely self contained inertial system, 
would provide world wide coverage, all weather ooeration, 
an unlimited number of users, availability of frequent 
fixes, immunity to interference, passive user ooeration, 
and high positional accuracy. Certain military operations 
reouire Positional accuracies of 0.1 nautical mile or better. 

Present long range navigation systems and their 

, [ 2-1 
characteristics are presented in Table 1.1, where long 

range is defined as having a position fixing capability at 

a range of approximately 400 nauticaL miles or greater. Of 

these, the inertial and acoustic doppler systems are 

essentially sophisticated dead reckoning systems with 

potential wo^ld coverage and position indicating capability. 

However, at the present time and state of the art, the use 

of a navigational satellite system appears as the best 

solution to the requirements outlined above. 

1.2 CONCEPTION OP A DOPPLER SATELLITE NAVIGATION ^ V S"EM. 
The concept of a doppler navigational satellite system was 
by-product of the observation of the Soviet Union's first 
artificial earth satellite. Sputnik I, in the autumn of 
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TAEL'S 


1.1 - LONG-’ 


2AFGS FAVTGATION 


S'VSTRIVS 


System 


Estimated 
Useful Range 


Estimated 

Accuracy 


Frequency 

Range 


Celestial 


World-wide 


3nm 


- 


Star 

Tracker 


World-wide 


0 • 3nm 


- 


Consol 


700-1200nxn 


0.3-0. 6° 


1 90-1 94k c 
257- 263k c 


RDF 


300-600nm 


2-6° 


200-41 5kc 
2R5-325kc 


Oecca 


200-300nm 


0 . 23-1 .Onm 


90-130kc 


Loran-A 


700-900nm 


0.5-3.0nm 


1 . 75-2mc 


Loran-C 


1200-1500nm 


0 . 2-0 . 5nm 


90-110kc 


Omega 


5000nm 


0 . 5-1 . Onm 


10-1 5k c 


Satellite 


World-wide 


0.1-2. Onm 


100-400mc 


Inertial 


3-12 hours 


2.0nm/hr 


- 


Acoustic 

Dopoler 


300 ft 
water deoth 


ltf of 0.2-lmc 

velocity 
l°/ 0 o f distance 



traveled 
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TABLR 1.1 (Cent) 



*!yp t ern 




Celestial 


Inexpensive equipment, easy to use; 
limited to 2 fixes/d-y, affected by 
clouds and bad weather which also 
limits fix availability to 50“/. 


Star Tracker 


Same as above except equipment 
very expensive. 


Consol 


Requires only a narrow band receiver 
to establish LPP. 


RDF 


World-wide coverage. U.S. coastal 
beacons. Capable of coarse fix only. 


Decca 


Rosy to use, insufficient coverage, 
limited to about ‘ 3 4-0 mile range at 
night because of sky wave interference. 


Loran-A 


Easy to use, insufficient coverage, 
subject to sky wave interference. 
Less accuracy with sky wave. 


Loran-C 


Rasy to use, accurate, insufficient 
coverage. Less accuracy with sky 
wave. Expensive installation. 


Omega 


World-wide coverage. Presently in 
advanced development. Will use three 
f^eouencies to yield both coarse and 
fine grids. 


Satellite 


World-wide coverage. Less sophisticated 
systems for coiner cial use now under study, 


Inertial 


Self contained, must be checked with 
outside reference. Accuracy decreases 
in a few hours at rapid rate. Very 
exoensive, requires expert maintenance. 


Acoustic Doppler 


Sophisticated high accuracy systems 
utilize gyro-compass, velocity input, 
and automatic track-plotting equipment. 
Limited to shallow water. 
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1957. Dr. W. H. Guier anii Dr. G. 0. Wei ff onbach of the 



Applied Physics Laboratory of Johns Hopkins Hniversity 
listened to the signal generated by Sputnik and noted 
the doppler shift in the transmitted radio frequency. The 
shift is given quantitatively by: 

(1.1) Af = -(p/ c )f s 

where Af is the doppler shift from the satellite transmission 
frequency, p is the range rate, c is the velocity of 

electromagnetic propagation, and f s is the satellite 

[3] 

transmission frequency. Hence, the measuring of doppler 
shift is comparable to measuring the range rate. This is 
the basis for the development of position determination 
equations derived in Chapter 3. 

Guier and Weiffenbach deduced, and later demonstrated, 
that accurate measurement of the doppler shift pattern 
would permit the determination of a satellite orbit. After 
the first analysis and use of the doppler tracking technique. 
Dr. F. T. McClure of the Applied Physics Laboratory, r en soned 
that the technique could be inverted and that the doppler 
curve received from a satellite whose position was known 
could be used to determine the position of the receiver 
on the surface of the earth. It seemed only logical to 
utili/.e this deduction as a navigation system. 

1.3 THH OP THAT TONAL RH0 TT IHHMT Tm . - - The Navy is 

constantly seeking improvements in the fundamental accuracy 



of navigation and in the operational character of navigation 

v . 

'V s 

syst em£ * Realizing that the d op pi er technioue provides a 

powerful tool for simple but very accurate naviration, tne 

development r»f a highly accurate operational navigation 

system using artificial satellites was made on* of the T'avy’s 

first space objectives. The Chief of Naval Operations 

established an operational reouirement for a navigational 

satellite as follows: 

Develop a satellite system to provide 
accurate, all-weather, world-wide navigation 
for naval surface ships, aircraft, and 
submarines . [d] 

This reouirement was generated by the need for a highly 
accurate position determination required by modern weapons 
systems. The satellite system was not developed to replace 
any of the present systems of navigation. Instead, it was 
develooed to provide another navigation aid for determining 
positions of extremely high accuracy anywhere on the surface 
of the earth. 

Accordingly, the Navy proposed to the Advanced Research 
Projects A-gency (ARPA) that a project be undertaken to 
exploit the doppler technioue for navigational purposes. 

The program began in 1958 as a feasibility study under 
contract with the Applied Physics Laboratory. This program 
was designated ''Project Transit 1 '. 
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2. DEVELOP VENT OF AN OPERATIONAL SYSTEM 



2.1 GENERAL CONFIGURATION AND OPERATION - - After 
the initial conception it did not take long to determine 
the configuration of an operational donpler satellite 
navigation system. The ba^ic system consists of the 
following five grouos of equipment: 

1. Satellites. 

2. Tracking stations. 

3. Computing center. 

4. Injection station. 

•5. Navigational equipment. 

The system is designed to employ four satellites in 
polar or near polar orbits at an altitude of 600 miles and 
a period of 110 minutes. With orbital planes separated by 
45 degrees a fix is available at least every 110 minutes. 

If the satellite oasses within 10 degrees of the zenith, 
the interval is extended to 220 minutes since longitude 
can not be accurately determined in this situation. At 
higher latitudes the freouency of fix availability increases 
considerably. Typical satellite orbits for an operational 

r5 l 

system are illustrated in Figure 2.1. The satellite 
transmits two rigorously coherent freouencies controlled 
from the same oscillator. This allows a first order 
correction to be made for ionospheric refraction. In 
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FIGURE 2.1 - TYPICAL OPERATIONAL SATELLITE ORBITS 



8 



addition the satellite transmits orbital information am 
time signals synchronized with 1TT-2 C . The time signals 
are obtained from a clock which is regulated by the 
satellite oscillator. 

The overall plan of a typical system operation is 

L4.1 

illustrated schematically in Figure 2.2. t navigational 
satellite is shown at three different time periods during 
an orbit. During time period T-j the satellite signals are 
received by a ground tracking station as the satellite 
passes within radio line of sight. As each tracking station 
receives the transmitted signal, it determines the apparent 
doppler shift and reads out the satellite memory and time 
marker. Those data are then transmitted in digital form by 
teletype to the computer center. 

After receiving a sufficient number of pass reports, 
the comouter center uses the doppler data to compute an 
improved satellite orbit and then oredicts the satellite 
oosition for every two minutes in the nexu 16 hour period. 
When this is conroleted, the computer generates the parameter^ 
of the orbit from which the position of the satellite at 
any point in time can be determined. The computer also 
analyzes the time errors to give the clock rate, and 
determines the correction to both the clock setting and the 
clock rate. All of this information is chen sent to the 
injection station. 



9 



Project TRANSIT 

TIME T„ 







J 



SCHEMATIC SYSTE/A OPERATION 

FI6URL ?. .2, V 



Dor in? time period Tp, when the satellite is within 
range of the injection station, the satellite memory is 
erased, the new data is inserted into the memory, and the 
clock is reset and regulated. The satellite immediately 
retransmits this message to the injection station where it 
is compared with the information originally transmitted. If 
the readback indicates that the injection was unsuccessful , 
the message is automatically retransmitted to the satellite. 
When accurate storage in the satellite memory is verified, 
the memory is locked by a time clock. Additional information 
is rejected until approximately twelve hours later when 
the satellite is again in optimum range of the injection 
station. 

After receiving the message from the injection station, 
the satellite continuously transmits the parameters of its 
orbit and its predicted orbital position at two minute 
intervals. During time period the navigation eouipment 
aboard a ship receives the orbital oarameters from the 
satellite and records the dotmler shift for at least three 
2 minute time intervals. Comnuting equipment uses this 
information to calculate the ship's true latitude, longitude 
and time. 

2.2 FEASIBILITY DETERMINATION. - - The first phase 
of the Transit Program was orimarily concerned with proving 
the feasibility of such a dopnler satellite navigation 
system. This phase was under the direction of ARPA. 
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During this period four experimental satellites were 
designed. They were designated the 1A, IP, 2A, and 2B 
respectively. Seven ground tracking stations were nut 
into operation and others were in the process of being 
designed and constructed. An experimental computation 
center was placed into operation and a time and freouency 
standards laboratory was established at the Applied Physics 
Laboratory. The requirements of the injection station 
were determined and the design of the station was undertaken. 
Extensive mathematical analyses of computational procedures, 
ionospheric refraction, geodesy, orbit determination, and 
position determination were initiated. 

On 17 September 1959 the Transit 1A satellite was 

launched by means of a Thor-Able vehicle. The satellite 

A'\ 

in it's flight shroud is shown in Figure 2.3. Although 
the satellite failed to obtain orbit, sufficient data were 
obtained to verify the feasibility of satellite tracking 
and navigation by doppler analysis. However, objectives 
concerned with the verification of the satellite design 
were not achieved. 

The Transit IB satellite was similar to the 1A in 
appearance and function. The satellite consists of a 
shell divided into two hemispheres 36 inches in diameter, 
a central support, and an instrument troy. The shell is 
a lamination of fiberglas with a honeycomb plastic filler. 
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This provides a strong, highly insulated, nonmetallic 
structure. White paint on the shell surface hids in 
temperature control. The central support, also made 
of laminated fiberglas, connects the hemispheres of the 
shell and supports the aluminum instrument tray. The 
tray is tied to the shell by nylon lacing which absorbs 

the shocks produced during launching. A cutaway view 

r. i 

of the Transit 13 satellite is given in Figure 2.4. 

The tray contains two transmitting systems, silver- 
zinc and nickel-cadmium batteries, two command receivers, 
and a de-spin system. A radiation and insulation shield 
is placed above and below the tray. Two banks of solar' 
cells used for charging the nickel-cadmium batteries and 
two de-spin weights are mounted outside and around the 
eouator of the shell. Broadband antennas in a spiral 
pattern are painted on the outside of the shell. 

The tray also contains an infrared scanner which 
was designed for snecial studies by the Naval Ordnance 
Test Center. A spin rate of approximately 180 rom was 
imparted to the satellite for this test. Aft er a ore-set 
time the de-spin system stops the spinning of the satellite 
since spinning renders doppler data useless. 

The Transit IP satellite was successfully launched 
on 13 April 1960 by means of a two stage Thor- Able-P tar 
vehicle. The elliptical orbit has an inclination of 51.3 
degrees, a perigee of 235 miles, and an apogee of 475 miles 
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The satellite period is 95.5 i inutes. although the 
satellite was not injected into the planned circular, 

500 mile orbit, the ground tracking end computation oroved 
completely successful. Analysis of telemetered data 
indicated that the satellite components functioned 
satisfactorily. All electrical and mechanical aspects of 
the satellite were proven including the solar cell operation, 
command off -on operation, de-soin devices, and temperature 
stability. Navigational experiments were conducted which 
gave position determinations to within - 1 - mile of first order 
geodetic control points. All calculations used in orbit 
and navigational fix determination employed the first order 
ionospheric refraction correction obtained from comparing 
the two coherent transmitted freouencies. The accuracy 
obtained indicated the validity of this two freouency 
techniaue. Investigations of the earth's gravitational field 
and geodetic surface were also conducted. 

The Transit IB satellite ceased radiating on 11 July 
1960. However, the results of this experiment achieved 
the goal of the first phase of the Transit research program 
by conclusively proving the system feasibility while 
demonstrating the ability to perform dopoler satellite 
navigation with an accuracy required for most military 
applications . 

2.3 PHOTOTYPE DYV dLOPM YNT . - - With the success 

of the Transit IB satellite the Transit Program entered 
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the second ohase of development in early May 1960 when 
the management of the Program was formally transferred 
from ARP A to the Bureau of Naval V/e aeons . This phase 
consisted of the development of an engineering prototype 
of the operational satellite navigation system. Cense uently, 
a series of research and development satellites were 
constructed and placed into orbit. The basic characteristics 

and orbital information of these satellites is given in 

[s] 

Table 2.1. effort was continued in improvement of system 
design, simplification, and reliability. Comprehensive 
testing in environmental chambers, and for vibration, shock, 
and acceleration were conducted on the satellite and it's 
components. -A complete data storage and readout system was 
developed including the design and construction of an 
injection station. Tracking procedures and equipment 
were improved along with the development of shipboard 
navigation equipment. Continued research was conducted in 
the areas of gravity, refraction, and geodesy resulting in 
improved accuracy. 

The Transit 2A satellite was placed into orbit on 
21 June 1960. It differed from the 1A and IB in that it had 
more solar cells, an imoroved telemetering system, and an 
electronic clock. The clock, running from the satellite 
oscillator, orovided the accurate time correlation between 
the tracking stations which is required for accurate 
determination of satellite orbits. 
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TABLT? 2.1 - TRANSIT AKJ) DKVU.OP 

SAT?XLITBS LAT T HOH"vD STRING 1960 AND 



M^NT 

1961 






Nov IB 


Nov 2A 


Nav 3B 


Nav 4A 


Nav 43 


TRAAC 


LAUNCH DATE: 


APRIL 13. 1960* 


JUNE 21. 1260 


FEB 21, 1961** 


JUNE 29. 1961 


NOV 15. 1961 


NOV 15. 1961* ** 


* 8" SYSTEM; 
FREQUENCIES "C SYSTEM; 

*T SYSTEM; 


162-2l6m< 

54-324m< 


162*2l6mc 

54*324mc 


162-216mc 

54-324m( 


54-324mt 

lS0-400mc 


54*324m< 

l50-400mc 


54-3 24 m( 


RMS FREQUENCY NOISE AS *‘B" SYS. 162-216mc 
MEASURED AT TRACKING "C" SYS. 54 324m( 
STATIONS; "V SYS. 150-400mc 


9 

5 Port* In 10. « 
5 Ports In 10 


9 

1 Port In 10 .q 
7 Port* In 10 


5 Paris In 10 lg 
5 Ports In 10 


10 

2 Parts In 10ig 
2 Ports In 10 






MEMORY; TYPE 

; CAPACITY 


Non* 


None 


Mognelic Shift 
Regis) er 
3B4 Bils 


Delay Line 
2049 Bils 


Ferrile Cave 
1344 6.1s 


None 


WEIGHT 


265 lbs. ; 


| 223.3 lbs. 


290.3 lbs. 


175 1 lbs. 


198 lbs. 


233 lbs. 


INITIAL PERIGEE ALTITUOE: 


37B km 1 


| 621 lm 


178 km 


383 U 


950 km 




INITIAL APOGEE ALTITUDE; | 


754 km 


| 1070 U 


97B km 


994 km 


1111 km 




INCLINATION; 


51.3 deg. 


66.7 deg. 


2B.4 deg. 


66.8 deg. 


314 deg. 


32.4 de?. 



* Note: Xav IB ceased radiating on July 11, 1960. 
** Nav 3B entered the atmosphere on March 30, 1961. 
*** Launched pickaback on Transit 4B. 
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The Transit 3B entail it a, launched on 21 February 1961, 
incorporated an additional electronic unit in the form of a 
complete data storage system capable of storing a small 
number of bits of digital information . A storage loading 
techniooe and readout system was used to insure correct 
data insertion. This provided the first flight test of 
the data system including the injection station. 

T "nfortunat ely, the 3B satellite only remained in orbit a 
little over five weeks and the experimental results were 
quite limited. 

The Transit 4 A and 4E satellites wore launched on 

29 June 1961 and 15 November 1961 respectively. A cutaway 

r cr- 

view of the Transit 4A satellite is given in Figure 2.5. 

The 4 series was extensively engineered to achieve an 
extremely stable oscillator, small size, and increased 
reliability. Crystal oscillntoi’s of high accuracy were 
carefully aged and selected. Their operational temperature 
was actively controlled by a small heating coil contained 
in the thermal inertial element. This emphasis on constant 
operating temperature yielded oscillator stabilities 
up to one part in 10 10 for short intervals. Such accuracy 
gives the satellite system a capability for extremely 
precise position determination. 

The internal equipment was miniaturized to reduce 
the satellite size and weight. This was to enable the 
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use of tin inexpensive Scout launch vehicle. The outer 
shell of solar cells is collapsable and folds back against 
the satellite outer shell assembly. This greauly reduces 
the size of the satellite as shown in Figure 2.5. Once 
the satellite was out into orbit the outer shell of solar 
cells was to be mechanically extended. In practice this 
was never achieved since the launch vehicle was always large 
enough to allow the satellite to be fully extended within 
it's shroud. The A series also contained an experimental 
radio isotope (FIPS) power source. 

In order to increase the satellite reliability, plated 
circuits and connection welding techniques were used through 
out. In addition, entire subsystems were imbedded in plasti 
to introduce the necessary ruggedness, strength, and 
temperature inertia. 

All the experimental satellites- so far have had to 
use isotropic antennas because they were not stabilized 
with respect to the earth. Since no antenna is truly 
isotropic the signal in some directions is much weaker 
than in others. As the satellite rotates, it sometimes 
sweeps a direction of low signal strength across the 
ground receiver and the signal lock can not be maintained. 
This accounts for almost all o e the satellite passages 
that can not be used. 

The Transit Research and attitude Control (TRAAC) 
satellite was launched pickaback by the same vehicle as 
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t ho 4R satellite. Its primary objective was testing a 
method of gravity attitude stabilize tion. It uses the 
gr .dient of the gravitational field acros? the dimensions 
of the satellite in order to keep one face of the satellite 
always toward the earth. This allows the use of a directional 
antenna and avoids the problem of low signal strength in 
the antenna pattern. 

As a result of the data gathered from the research 
and development satellites the basic characteristics and 
performance of a satellite navigation system were validated 
and a final prototype operational satellite was developed. 

The first of the operational Transit 5 series was launched 
in 196 2 . An artist’s conception of the Transit 5A satellite 

C ' 

is illustrated in Figure 2.6. Details pertaining to the 
operational satellite and system are not readily available 
due to security restrictions. However, this scstollite 
incorporates the design features proven in the previous 
experimental satellites. Conspicious additions include 
sola v> cell paddles and a damping system. The solar cell 
'paddles can be angularly oriented with respect to the sun 
and thereby insure a nearly constant output of cower. The 
damoing system is comoosed of a long boom with a small 
weight attached to a spring at the end of the boom. The 
large displacements of the spring provides the necessary 
damping of angular oscillations reouired to keep the 
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FIGURE 2.6 



P R OT 0?vp R OP E° AT I DUAL 
TYPE 5 A 



SATELLITE 



satellite gravity stabilized. 

Since July 1962 the TT . S. Naval Observatory has 
maintained equipment for the recovery of time signals from 
Transit type satellites to study the problems associated 
with providing continuous synchronization of TT T-2 C . precise 

uime is required for highly accurate orbit and position 
determination. 

A later model of the prototype operational satellite 
was designed to use nuclear energy as a main power source. 
When this satellite was placed into orbit it was the first 
operational satellite to use this new energy in space. To 
date, a number of prototype satellites have been placed 
into orbit. 

2.4 THE OPERATIONAL SYSTEM. - - In July of 1964 
three of the prototype satellites were incorporated into 
an operational system. A complete operational ground network 
v/ as set up and is entirely separate from the research and 
development network. The research and development network 
is 'still in existance and is designed to provide basic 
information on the size and shape of the earth and it’s 
gravity field. This information is periodically processed 
through a computer at the Applied Physics Laboratory and 
occasionally new force terms are derived. The new terms 
are then used by the operational ^y^tem to improve orbital 
updating and predictions. In July 1965 the name of the 
Transit System was officially changed to the Navy Navigation 
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The use of Transit is no longer authorized 



Satellite System, 
within the Navy. A discussion of the ooerationnl system 
and it’s components is nresented in Chanter 4, following an 
investigation of the mathematical theory of satellite 
navigation presented in Chapter 3. 



3. NAVIGATIONAL SATELLITE THEORY 



3.1 GEYEEAL. - - We will assume an artificial 

satellite in a known polar or near polar orbit transmitting 

on two extremely stable coherent freouencies in the 100 to 

400 megacycle range with added uniformly spaced time signals 

At a receiving station an apoarent change in the satellite 

transmission freouency is observed.' This is caused by the 

motion of the satellite and the receiving station. The 

amount is prooortional to the velocity of approach or 

recession. The exact amount of this shift depends on the 

exact location of the receiving station with respect to the 

S 

path of the satellite. This is illustrated in Figure 3.1. 
The navigation equipment integrates the heterodyne frequency 
v/hich is in the tens of kilocycles range, between the 
satellite transmission and a reference oscillator in the 
navigation equioment . We note from Figure 3.1 that the 
sign of the frequency shift changes from plus to minus at 
the center of passage. ' The beat note does not distinguish 
the sign of the freouency difference. In order to get 
correct measurement «? it is required that the beat freouency 
does not pass through zero. Therefore, the local oscillator 
frequency is offset from the satellite oscillator by 80 part 
in 10^. This offset exceeds the maximum doppler shift 
that the beat note never changes sign. Uocn receipt of each 
time signal the number of cycles counted during the previous 
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FIGURE 3.1 



EFFECT OF OBSERVERS LONGITUDE DISPLACEMENT 
FROM SATELLITE SUBTRACK. 



time interval is recorded. The na viga t ion eauipment uses 
the higher of the two freouencies to determine count and the 
comonriaon between the two for a first order correction for 
ionospheric refraction. It is then assumed that the recorded 
cycle count is refraction free. 

3.2 REFRACTION C^RR'XTION. - - Accurate do-pler 
measurements have utilized the dispersive nature of the 
ionosohere to eliminate the refraction contribution which 

Zsl 

is expressed as a series expansion in the inverse frequency. 

f A ( t ) B(t) C ( t ) D ( t ) 

(3.1) Af(t,f )+ — /Oft) = - - ~ 

c 1 f s f| fj 4 

By receiving two coherent freouencies transmitted from the 
satellite the freouencies may be substituted into (3.1) and 
the first unknown coefficient 'A' determined at each time 
value t , where p is the range rate and c the velocity of 
electromagnetic propagation. Y/ith the use of only two 
freouencies the resulting cycle count still contains all 
contributions that depend on the higher powers of (l/f ). 

In operational systems two frequencies are generally 
adequate for most precision calculations. The tracking 
station receiver uses analog methods to automatically make 
the refraction correction before the dopnler data are 
digitized. A further discussion of ionospheric refraction 
is presented in Section 5.2. 
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The 



o.o DYT^RY IKnT IOK 0? TrfY I jTYGR.-.TiD CYCL3 C..UM. 

The satellite transmits a time signal at T^ and at T 2 « 

transmission times from the satellite to the receiver are 
t-|_ and tp. If we denote: 

f = satellite transmission ^reouercy. 

f r = reference frequency. 

f g = observed satellite transmission freauency 
affected by doopler shift. 

IC — integrated cycle count. 

we can derive the following relationships. 

The heterodyne frequency observed at the receiver is: 



(3.2) 



f h = f" - f 
n s r 



and the integrated cycle count is obtained from: 






(3.3) 



IC = 



(f‘; - f\.)dt 



T 



Vtl 



Solitting the integral and evaluating the last term: 






(3.4) 



( f " ) dt 
v s ' 



- f. 



( w 



(T 1 -f-t 1 ) 



Tn- 



The integral of (3.4) is nothing more than the number of 
cycles transmitted by the satellite during the uniform 
time interval T 2 - T-j_ . 
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we can say: 



Therefore, 



( 3 ♦ b ) 



r 



(I ~~ j 

ip' Op 



r 



2 



( f b dt 



m jl. -J- 

A l’ l 



i <V dt 

/ T n 



Substituting the result of (3.5) into (3.4): 



(3.6) IC - f g (T g 



V - f r <VV - < T l +t l> 



By rearranging terms (3.6) can be written: 



(3.7) 
Where : 



ic = Af AT - f A 



A f 

1 T 
At 



- T, 



— tn — t- 



Integration allows the use of range differences (see 
Souation 3.9) instead of range rate in computations. 

3.4 RFLnTION OF DISTANCE AND INT-2GTIAT3D CYCLR CHTFT. - 
The basic distance formula states that distance is equal to 
velocity times time. hence, the difference in transmission 
times from the satellite is proportional to the difference 
in range at and T g . Solving (3.7) for At: 

Af AT - TC 



(5.8) 



At - 
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If v/e lot r eoual the rsnpe to the satellite we can express 

(3.8) in terms of distance. 

c( A f \T - IG) 

(3.9) Ar — 



With this basic eouation relating distance and integrated 
cycle count we can derive equations Ar position 
determination in a rectangular coordinate system. 

3.5 EQUATIONS FOR POSITION T3TFRT/INATI0N . - - 



If we denote: 

r — true distance to the satellite. 



R — computed distance to the satellite from 
a DR position. 

X ,Y ,Z — assumed ship coordinates in a right 
n n hand orthogonal coordinate system. 

(obtained from some other navigation 
system) 

X?,Y?,Z S — satellite coordinates in the same 
system. 

i — 0,1,2, N - where N is eoual to the number 

of time ooints. 



Then : 

(3.10) r d.R 

Using assumed values for the variables as opposed to the 
true values in (3.9) and differentiating (3.9) and (3.10) 
with respect to the variables we arrive at: 



(3.11) A r - R, 



M. - R i+ d Vl - A= 7 ( “ f n AT - IC n> 

r 



c c 

— (ATdAf n ) dIC 



n 



■ r 
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Rxpres «inp the R snd dR terms in terms of our orthogonal 



coordinate aystom: 



(5.12) 



(3.13) 



R. 

l 





- X 

n 





ihx) 2 -| ( Y s - Y - iiy) 2 
i. i n 

iAz) 2 p 

X n - i.ix)dX n -; (Y* - Y n - iAy)dY n 



+ (Z* - Z n 



z) dZ 



n 



B'here Ax, Ay, and Az account for navigator motion between 
observations. Simi.liar enuations can be developed for R^ 
and dR^^ . By substitution of values obtained, above into 
(3.11) Ar can be expressed as: 

(X? - X n - iAx)dX n -|- 
'•i 



(3.14) Ar = R.^ , - R.r 
i+l l 



(Y. - Y n - iAy)dY n (zj - Z n - ilz)dZ^ 
1 

- — (X s - X - ( i-1 )Ax) dX n -j- 
R. - i^ 1 n 



i-1 



( Y . 



_ V 



i-rl 



n 



- (i+l)Ay)dY n n-(Z^ 1 - Z n - (i-il)Az)dZ n 



TTsing this exoression for Ar in (3.11) and rear -anging terms 
we can form an observation type equation where v is the 
residual formed by the adjusted value of Ar minus the 
observed value of Ar. 

= Ar a " Ar b 



v 
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If we denote the coefficients of (3.14) as follows: 



— coefficient 


of 


d.X n 


b.. — coefficient 


of 


dY n 


c^ zr coefficient 


of 


fiz n 



the eouation takes the form: 



c 



(3.15) 



v • 



a.dX 4 - b dy 4- c dZ AT dAf. — 

1 n i n ] *1 r> n 



(7 



i-i-l 





IC b ) 



■“r 

In this equation we assume our cycle count is errorless 
and therefore neglect the dIC term. The navigator may 
also assume that the frequency of his reference oscillator 
is unknown hut constant during the short interval of 
observation. In the equation this is treated as the unknown 

dAf. The above eouation contains four unknowns (dX_, dY , 

n> n* 

dZ n , dAf n ) • If we have four satellite observations we can 
solve for the unknowns and thus locate our position. 

Eauation (3.15) has been called an observation type 
equation instead of an observation equation because the 
residual is formed in terms of A r instead of 10. In 
Section 3.7 the actual observation eouation is derived 
where v is the residual formed by the adjusted value of IC 
minus the observed value of IC . It will also be shown that 
for navigational purposes Equation (3.15) may be used as the 
actual observation eouation in computations. 



If we assume the earth to be approximated by an ellipsoid 
of revolution then: 



( 3 . 16 ) 



r2 

'n 



O 

Y° 

n 



Z 



2 




V/e can reduce the number of unknowns to three by expressing 
dZ n in terms of dX n and dY^. Differentiating (3.16): 

' X n x n \ o' 

(3.17) dS - dx n - — dY n : l - 

Z„ Z I 

n n 1 i 

The unknowns are now dx n , dY n , and dAf. Now only three 
observations of a satellite are renuired to fix a position. 
By substitution of the value of dZ n from (3.17) into (3.14), 
multiplying and gathering terms, we can form an observation 
equation similar to (3.15). 



(3.18) 



v i 



R . , , - R.4- 
i+l i ‘ 



*i \ 



(X? - X - iAx) 



X 



n K 



(Zi - z n - iAz) (1-e^) ( ) 



Z n , R i+1 ; 



(Xirl " X 



n 



X \ 

( i-KL )Ax ) — j— ( 2 7 - Z n - ( i-rl )Lz) (1-e 2 ) ( --) 

Z n / 



dX 



n 



*i i 



(Yf 



Y n - iAy)-|-(Z| 



Y 

7, - iAz)(l-e 2 )(-— ) 

7 

"n j 
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1 f 
- — (Y? 



R i-tl 



Ul' y n' U+U^yJ + U^- Z n - (i+1)/ z) 



n 



y ! c c 

(i_ e 2)( — £) : dy (/ f At - ic.) — At dzf 

11 ~ D o 



z 



n 



fv 



f v 






This equation can be reduced to the form: 



(3.19) v. = a.dX -4-b.dY -f c.dZ.f-fL. 

i i n i n i n i 



n 



where a^ — coefficient of dX n 

b^ coefficient of dY n 

c^ ~ coefficient of dAf 
1 o n 

L i = <h +1 - V - — Af b iT - I0 b 

1 r 

In practice we constrain our position to be on a sohere 
since we are dealing with small quantities. 

Equation (3.16) becomes: 



(3.20) 



O O O O 

T J -f- Y -f- Z° = R 2 
n n n e 



and (3.17) becomes: 



(3.21) 



dZ n - - 



n 



Y, 



,e v 

-iV 



n 



n 



dY 



n 



“*n n 

The maximum difference in the dZ n term resulting from 
this approximation can be determined as follows: 



(3.22) 

and 



dZ 



n 



( ellipsoid ) 



- (1 - 9")dZ n 



( sphere ) 



d.Z 



n 



- dZ 



( ellipsoid) 



n 



- e 2 d Z 



( sohere) 



n 



( sphere ) 
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f 













Therefore, the difference in (3.17) due to constraint to 

O 

n sphere is -e a1Z n ( pp^ero) • The terra dZ n bn? it's greatest 
value at the eouator. Assuming that the navigator position 
is known to an accuracy of one nautical mile (-■- 1829 meters) 
and using the satellite derived value of e : ’ from the Kaula 
Ellipsoid of 1965, we find the maximum difference in dZ n 
to be ; 

AdZ n — (0.00669) (1829) or czl2.2 meters. 

This result is considered negligible for the purpose of 
navigation. By using the spherical approximation the 
(1 - 9^) terms in the coefficients a^ and b^ are eliminated 
and the observation equation retains the form of (3.19). 
This is the general satellite observation equation. 

3.6 REDACTION OF OBSERVATION E STATIONS. - - Three 
observations of a satellite will yield th^ee observation 
equations. In practice these are used directly in the 
determination of the three unknowns, ’"hen there are more 
observations than unknowns the method of least squares may 
be applied to find that system of values which is most 
probable from the observations themselves. Arbitrarily 
assigning a weight of unity to each observation and adding 
the square sums of the coefficients we obtain: 

(3.23) [vy] = jha] dX^ +. 2 [ a b] dX n dy n + 2 b<? dX n dAf n +2 [aL] dX n *h 
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[bbj dY^-g 2 |bcj dY n dAf n a |oL] dY n + jcc] dAf^ + 2 jctf d£f n -h{LL] 



where: [u§ = a^. -g a i+1 a. +1 

By the method of least souares the most likely solution is 
when the sum of the squares of the residuals [v\f is a 
minimum. To obtain this condition we set the partial 
derivatives of (3.23) with respect to the three variables 
equal to zero. This results in the formation of three 
normal equations. 

[aa)dX n +l?xb]dY n Biacn dAf n +gg - 0 

(3.24) &b]dX n +-!bbJ dY n + [be] dAf n + bLj = 0 

jac; dx n + ibc: dY^ [ccj d£f n + jcL| - 0 

Equations (3.24) are normal equations reduced from satellite 
observation equations. Regardless of the number of 
observation equations we can always reduce them to a number 
of normal equations equal to the number of unknowns. The 
solution of the normal equations yields unique results for 
the unknowns dX n , dV n , and dAf • Taking the values for dX n 
and dY n we can solve for dZ p . Addin# these corrections to 
the original values for X^, Y n , and Z n gives the rectangular 
coordinates of our adjusted position. 

(3.25) X a = X n + dX n 
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= r n + iy n 



Z = Z + <iz 
a n n 



The value obtained for dAf n may be used to to update the 
freouency of the reference oscillator in the navigation 
eouiprrient . 

It should be again noted that the use of the least 
squares solution is advoided in practice by using only 
three observations to solve directly for the unknowns. If 
more than three observations are obtained the excess 
observations may be neglected or they may be combined to 
oroduce a total number of three. This is the most practical 
solution since the advantages of the least squares solution 
are minimized by the small number of observations that are 
available during one satellite n.ass. During a maximum pass 
time of approximately 1ft minutes only seven observations 
can be made. at the oresent time a decrease in the time 
interval to gain more observations is not practical due to 
comoutational time requirements and restrictions of the 
satellite, memory and transmission system. 



5.7 DERIVATION OF THE OBSERVATION EQUATION • - - In 
Section 5.5 an observation type equation was formed where 
v is the residual formed by the adjusted value of Ar minus 
the observed value of Ar. In the theoretical case the 
observation eoua tion ■ should be derived in terms of the 
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observed Quantity. Such an observation equation will now be 
develooed. in terms of the integrated cycle count where the 
residual is formed by the adjusted IC minus the observed IC . 
It will then be shown that for the purpose of navigation 
the actual observation equation can be replaced by Equation 
( 5 • 1 o ) « 

TJsing the notation of Section 5.5 and denoting: 

IC — adjusted integrated cycle count. 

IC^ r: observed integrated cycle count, 
we can exnress IC as a function of the unknown parameters. 

(3.26) IG a — f(X a , Y q , Z a , Af n ) 

and : 

(3.2V) IC a “ IC b -|- v or v — IC g - IC b 

The adjusted parameters can be expresses as: 

X a = X n+ * 

Y a = Y n + y 
Z a - Z n 1 z 

Af a = Af n+ p 

where x, y, z, and F are corrections to the assumed, navigator 
parameters. Assuming that the corrections are small comoared 
to the parameters we can exoand (3.27) using a Taylor’s 
series whore terms of the second order and greater are 
negl ected . 
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I 






, <bf >f 

(3.28) IC, -f V - f(X , Y , Z , A f ) -j- ~ X .J. y -'r z 



n n n 



n \y v Y '■ 7 

0X n - Y n < Z n 



).f 



c Af 



P 4 



n 



ivnsre : 



f(x n , y , z , A f ) - ic, 

n n n n n 

The form of the actual observation equation for doppler 
satellite observations can be expressed as: 



(3.29) v . 



0 f 

^x. 



of 

x \ — y -t 

> y 

o J 



\ f 

CJ , 
z -V- 



\ f 
c* x 



— F 4- IC - IC K 

c') Z oh.f b 

n n 



n n 

In order to evaluate the coefficients of the corrections we 



note from (3.9) : 



IG =ATAf -Ar 

c 



and by substitution: 

(3.30) IG = AT/\f n 






- x n> 2 + (Y i + l - Y n> 2 



» 1 



+ < Z Ll - Z n> 



,(x? - x n )2 (y = - y n ) 



, fl 



< Z I - *»> 



\ 2 



The coefficients 'can now be determined from (3.30). 



0 ic 



c>x 



n 



! i 



i J .-l 



] / 1 

(X i + l - X n)+ — (A " V 



1 



J 



dX 



n 
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1 



1 



( Y? - Y ) j dY 
' 1 n' ~n 



^IC 

* Y n 



J 7 ii-l 



< Y Ll 



Y n>, 



i \ Hj 



^ IC 

d Z n 



( 1 



T~ U i+l - 7 -n>, 

A-KL 



Ui (Zi " Zn) '_ 



dZ 



n 



^ IC 

At d Af 

*• Af 



n 



If we denote: 



aV r= coefficient of dX n 



coefficient of dY 



n 



c m coefficient of dZ n 



we can write (3.29) in the form; 



(3.31) v. = a'JdX n ^b':'dY n -hc*dZ n -t- ATd Af n d- IG n - IC b 

c 

If we multiply (3.31) by a constant equal to we 

f 

Vk 

arrive at : 



(3.32) 



v i = V 1X n+ b i d V- c i 



dZ 



n 



ATd Af 



n 



(ic n - I0 b ) 



The coefficients of this equation are equal to the 
coefficients of (3.15). Prom (3.15) we obtain; 



Ar n n (R i+1 - K ± ) ~ 



(A T Af - IC n ) 
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and : 



Ar b= -j- (AT/_f b - ic b ) 

*r 

(3.15) can now be expressed as: 

c c 

(3.33) v i= a i dX n -V-b J ,dY n + c.dZ n ATdAf n - (IC n - IC^) 

H At (Af - Af ) 

f* Tl O 

r 

From Equations (3.32) and (3.33) the relation between the 
actual observation equation and the observation type equation 
based on terms of A 11 can be expressed as: 

C Tp A y» C 

(3.34) v = v r -t- — At (Af - Af, ) 

& i ir* n o 

r r 

and : 

(3.35) v f r = K v i°+ KAT (A f n - Af b ) 

c 

where K is a constant equal to . 

f 

r 

Feuation (3.35) shows that vA r is related to v^"' by a 
constant except for the second term on the right side. It 
will now be shown that this second term can be neglected 
for the purpose of navigation. From (3.6) we note: 

Af - A - f r 
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and : 



Zlr 



n 



f 

s 



1 

r n 




From this it can be determined that the term (/\f“ n - Af^) 
is eoual to (f , - f ). V/e will now determine the 

aooroximate maximum magnitude of the second term. Knowing 
that f has a drift rate of less than 2 parts in 10^® per 
day and assuming that we use the solution of dAf n from 
previous observations to correct f , we can keen f r to 
approximately 2 parts in lO^'. Therefore, the term 
( A f n - A f^) has an approximate maximum value of (- 
or ±2.8 parts in lo"^A Using a value of 186,000 miles/sec 
for c and 400 mc/sec for f the approximate maximum value 
of the second term of (3.35) can be determined. 

186,000 2.8 

kAT( Af - Af.) hh (120) r-7- (4 x 10 b ) 

n b 4 x 10° 10 1U 



or C±i0.0062 miles or ±12 meters. 



Since this value can be considered negligible for 
navigational purposes v/e can write: 



( 3 . 36 ) 





Therefore, if v/e set 




IC~! 
r i J 



to a minimum it follows that 
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v^-' v> v^ r also goes. to a minimum. Hence we are able to 
use (5.15) in place of (5.31) in our eomoutst ions . 

3.8 EARTH RECTANGULAR COO^UTUAT IS . - - The adjusted 

coordinate values in a rectangular system are of little 
oractical use to the navigatior. Positional coordinates 
are desired in terms of latitude and longitude. While 
the calculations for coordinate determination and 
transformation are generally well known, they are included 
to provide a comolete development of a position determination 
from doopler satellite observations. 

In calculating the positions of navigational fixes on 
the surface of the earth it is necessary to consider these 
points as lying uoon some mathematical surface. The actual 
shape of the earth's surface is quite irregular and is not 
adapted to the purpose of coirmutation. For this reason it 
is necessary to adopt an assumed figure which will not d >Dart 
from the true siarface by an amount sufficient to produce 
serious errors in our results. The figure generally adopted 
is an ellipsoid of revolution which is defined by a given 
flattening (f) and a given semi-major axis (a). 

The position derivation given in the preceeding sections 
is with respect to an arbitrary right hand orthogonal 
coordinate system. We will now define this system with 
resoect to the ellipsoid of revolution assumed to represent 
the figure of the eai'th. The origin is at the center of the 
ellipsoid which is also the assumed center of mass. The 
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Z axis coincides with the mean rotational axis and the X 
axis intersects the Greenwich T.'eridian at the equator. The 
Y axis completes the right hand system. The coordinate 
system thus defined is illustrated in Figure 3.3. 

The relations between geodetic latitude and longitude 
and the defined XYZ coordinate system, which can be derived 
from Figure 3.3, are given in the following .formulas. 



cos 



X = 



(p CCS /\ 



(3.37) 



Y — 



a cos (p sin k )\ , 



W 



a(l - e 2 ) sin Qp 



’There ’7 is defined as : 

’.V — (1 - e^ sin""' Q|) ) 2 

3.9 DFTFRUINaTICN OF LATITUDE AND LONGITUDY. - - 

Using values Tor X,, Y , and Z from (3.25) and the relations 

a a a 

given in (3.37) we are able to determine the latitude and 
longitude on the surface of the given reference ellipsoid. 

The value for longitude can be solved directly from (3.37). 

(3.38) X — Tan - ^ 

X 

a 
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a — semi major axis 
b — semi minor axis 
9 s first eccentricity 
Cp = geodetic latitude 
= geodetic longitude 



FIGURE 3.3 - liARTH RECTAU^TLAR COORDINATE SYSTEM 
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The solution for* latitude is a little more involved. From 



Figure 3.4 we observe that; 

N sinCP 

(3.39) tan(T) = — - — — 7 

A (X 2 Y 2 )a 



Where N is the radius of curvature of the surface of the 
ellipsoid in a plane through the normal and at right angles 
to the meridian. This is also referred to as the radius 
of curvature in the prime vertical. A problem arises that 
N is a function of Q) . (N is equal to a/vV) To solve 

(3.39) we must eliminate (p from the right side of the 
equation. Going back to (3.37) we can express Z in the 
following manner: 

(3.40) Z = N sin (p • N e 2 sin(p 

and N sin Q = Z -\- N e2 sin(p 

Substituting (3.40) into (3.39): 

Zd-N e 2 sin(p Z N e 2 sin<p 

(3.41) tan(p — ~ 3. ~ ~ "i (14 ) 

(X 2 Y 2 ) 2 (X^ Y 2 )s Z 



Substituting N ( 1 - e 2 ) sinQp for Z in the denominator 
of (3.41) we arrive at: 



(3.42) 



<p = 



Tan ^ ( 



(x; 



Y 2 ) 

a 



1) 



1 - e‘ 



The values of latitude and longitude obtained from (3.38) 
and (3.42) are used by the navigator to fix his position 
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z 




FIGURE 3.4 - MERIDIAN SECTION IN EARTH COORDINATE 
SYSTEM. 



49 



on the surface of the ellipsoid. 

Reversing this procedure the navigator c^n easily 
calculate his position in earth r act angul ar coordinates 
from values of latitude and longitude by applying (3.37). 

3.10 INERTIAL COORDINATE SYSTEM . - - The derivations 

to this point have assumed that the navigator is given 
satellite coordinates in the defined earth rectangular 
coordinate system. In practice the navigator must calculate 
these coordinates from given orbital parameters which give 
satellite coordinates in an inertial system. The 
determination of these parameters is quite complex and is 
omitted from this thesis. Readers desiring this information 
are referred to Reference 7. 

The defined inertial coordinate system is illustrated 
in Figure 3.5. The origin of this system is the same as 
that of the earth rectangular coordinate system defined in 
Section 3.8. The Zj axis coincides with the pole of the 
orbital plane. The Xj axis passes through the perigee 
and the Yj axis completes the right hand system. 

3.11 SATELLITE COORDINATE DETERMINATION. - - All 
navigational satellite systems require accurate orbital 
information so that the satellite coordinates at any point 
in time may be calculated. Since it is not possible to 
predict the position of a satellite with respect to time 

to the required accuracy for more than a day it is impossible 
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FIGURE 5.5 - INERTIAL COORDINATE SYSTEM 
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To overcome this drawback the satellite contains an ephemeris 
table stored in a magnetic memory which is periodically 
updated with new parameters from the ground. 

The memory is divided into two parts. The first or 
fixed part contains the oarameters required to define the 
best elliotical approximation to the actual orbit. These 
parameters change slowly with time and are taken as constant 
over a twelve to sixteen hour period. The d escr eoancies 
between this orbit and the actual orbit are small. The 
second or ephemeral pact of the memory contains the small 
corrections to the basic orbit for each time point. The 
orbital parameters stored by the satellite memory are given 



in Table 3.1. From this data the navigator is able to compute 

the position of the satellite in inertial rectangular 

coordinates at the beginning and end of each time interval. 

To calculate the satellite coordinates , Y s , and Z s 

I I I 

the navigator goes sequentially through the following 

C3J 

calculations where t m is the mean satellite pass time. 



(3.43) 




t 



P 



epp = (pp^°V^p^ tm 



Ax — Am ^ ° V 



Mm = 3f0D (h- A mH-(J 0 £t 

360 



m 
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TABLF 3.1 



SATFLLI'TF CAPITAL P 



Fixed paramet 



t 



P 



I ; 
vj o 




W P * 




s i * 

Am * 




ors : 

time of first perigee after 0 n or 12^T T T. 
mean angular motion * 2Tf/poriod. 
argument of perigee at t . 
precession rate of perigee, 
eccentricity of orbit, 
mean semi-ma.ior axis. 

right ascension of ascending node at t 

tr 

precession rate of ascending node. 

cos { ■* [=■ inclination of orbital plane. 

sing 

change of mean anomaly for 1 hour, 
change of mean anomaly for 2 minutes. 
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ABLE 3.1 (Co nfc) 



phemoral parameters: 

t^ * time after integral half hour , A^. 

£ * correction to true anomaly at t^. 

£ *- correction to semi-major axis at t^. 
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M i = M m + (R - V S I’/a 

6"i = £> in Iv! i 

E i — Mi + 6"i +- cJ E i 
A i - A o + $ A i 

Then the satellite coorrlinates in the inertial system at 
time t^ are priven by: 

Xj - A j _( cos % - £_) 

(3.44) = Ai [sin 3 j _(l - £_ 2 )i] 

— 0 

These are the satellite coordinates in the defined inertial 
coordinate system at the i th time point. In order to be 
used for position determination in (3.18) the inertial co- 
ordinates must be transformed to the defined earth 
rectangular coordinate system. 

3.12 TRANSFORMATION OF INFRTIAL COOROTNATFS. - - The 
transformation of the inertial satellite coordinates to the 
earth rectangular coordinate system can be considered in 
two distinct steos. The first step involves two rotations 
to bring the inertial X^ and axes into the eouatorial 
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plane and the axis to coincide with the rotational axis. 

i t . 

The results of these rotations define the X T , Y , and Z 

I I I 

axes illustrated in Figure 3.5. The second step involves 
bringing the X-J- axis to coincide with the X axis. This is 
complicated by the earth's rotation about its axis. Since 
the earth rectangular coordinate system rotates with respect 
to the inertial coordinate system the second step takes the 
form of a rotation of the inertial crime system about the 
2 axis through an angle of - _/\_g - GAST . Upon completion 
of these steps the satellite coordinates will be in the 
earth rectangular coordinate system used in the position 
derivation. The matrix form of the required rotations is 
given below: 



Step 1 - rotation about the Zj axis through an angle 

„(°) 

of - followed by- a rotation about the 

Xj axis through an angle of - ( 



(3. .51) 

•" t" 



Y = 

I 



l Zi j 



0 



0 



0 cos ( - ( ) sin( - c ) 



0 -sin( -( ) cos( -( ) 



(o) 



(o) 



0 



0 



% - / \ ~ f r* 

cos(- . o ) sin(- ./ ) 0 X 



i I 



( O ) ( O ) 

-sin( - 4 cos ( • ) 0 Y- 



1 Z 

J. 
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Step 2 - rotation about the Z axis through an angle 

o f {- j\. N - gast). 

(3.32) 

cos(A h - GAST) -sin(A N - GAST) 
sin(A N - GAST) cos(A N ' GAST) 

0 0 

3.13 SUMMARY - - The theoretical development of 
eauations for position determination given in the proceeding 
sections illustrates how a navigator is able to fix his 
position on the surface of the ellipsoid used as an assumed 
figure of the earth from observations of a navigational 
satellite. Initial assumed coordinates X_, , v , and Z 
are reoxiired from some independent source. In addition, 
the navip-ator must determine the integrated cycle count 
from at least three satellite observations and be able to 
calculate the coordinates of the satellite at the time of 
these observations. The solution involves the determination 
of corrections dX n , dY n , and dZ n which are applied to the 
assumed initial coordinates. This gives the new position 
coordinates in a rectangular system. These coordinates 
are then expressed in terms of geodetic latitude and 
longit\ide . 
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4. OPERATIONAL SYSTEM COMPONENTS 



4.1 GENERAL. - - The Navy Navigation Satellite 
System has been operational since July 1964. It consists 
of three satellites, four tracking stations, an operations/ 
computer center, two injection facilities, and numerous 
navigational receiving units. The ground network is based 
entirely within the United States and is operated by the 
U. S. Navy Astronautics Group. The locations of the 

M 

components of the ground network are shown in Figure 4.1. 

The task of the Astronautics Group is to determine, 
process, and load orbital and time information into the 
memories of the operational satellites. These operations 
must be carried out within critical time limits. Each 
satellite is tracked every time it comes within range of 
one of the four tracking stations. The data received 
from the satellite is then sent to the computer which 
updates the orbit, predicts the satellite positions for 
the next 16 hours, and generates the orbital parameters 
in a coded message. This message must be sent to the 
injection facility in sufficient time to be injected into 
the satellite before it has used up all of the information 
previously stored in its memory. Each cycle requires 
about 12 hours per satellite. Since the satellite memory 
is limited to 16 hours of information, operations must 
meet an exacting schedule. 
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V/ahiawa, Hawaii 
Tracking Station 



4.2 THE SATELLITE. - - As was previously mentioned, 
details concerning the operational satellite are not readily 
available d\ae to security restrictions. However,' it can be 
assumed that the ooerational satellite employs the basic 
characteristics develooed by the Transit research and 
develoDment series described in Chanter 2. 

4.3 THE TRACKING STATION. - - The four ooerational 
tracking stations are located in Maine, Minnesota, California, 
and Hawaii. The tracking station equipment consists of a 
Westinghouse AN/BRN-3 type receiver, a doppler cycle counter, 
an ultra-stable oscillator, and communications equipment. 

The data obtained by the tracking station is put into standard 
teletype format consisting of two parts. The first part is 
called the header and contains the identifying information 
required by the computer. The second part contains the 
actual data. The data consists of a series of 11-digit 
data points. Each point represents one doppler reading and 
the time that the reading was taken. These readings are 
usually taken every four seconds. Generally 100 or more 
data points are obtained during the puss of one satellite 
within radio line-of-sight of the tracking station. In 
messages sent to the operations/computer center, the first 
five digits of a data point represent the time of day in 
seconds (UT) and the last six digits represent the number 
of microseconds it takes to receive a predetermined number 



GO 



of cycles from the satellite. This preselected value is 
such that the time approaches but never equals one second. 

In order to send the great quantity of data gathered 
by the tracking stations to the ooerat ions/computer center, 
high speed communication eauioment is required. A rapid 
data-link called TRADAT, which has a transmission rate 
five times faster than teletype, is presently being used. 
TRADAT is essentially error free because it incorporates 
self-checking circuitry. 

4.4 OPTRATIOKS/COlviPlITRR CTKTRR. - - The operations/ 
computer center is located at the headquarters of the 
U. S. Navy Astronautics Group at Point H'ugu, California. 

All data from the tracking stations and the computer pass 
through the operations center. This center at all times 
maintains the ooerational system status and the positions 
of the satellites os well as related systems information. 

It is the operations center's responsibility to insure 
that the tracking stations track each satellite and that 
sufficient data is provided to the computer to allow 
accurate orbit determinations and predictions. It schedules 
the computer center runs and insures that the orbital and 
time parameters are sent to the injection facility in time 
to meet the injection pass of a satellite. 

Four satellites in orbit will prodv.ce an average of 
84 messages per day from the tracking stations. The data 
in these messages pass through the operations center and 
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are sent to one computer center whose primary job is 
orbit updating and position prediction. The computer 
center has an IBM 7094 II/1410 configuration with a 
360/40 for command and control functions. The doppler 
data are fed into the computer which is programmed for 
orbital analysis. The machine computes a theoretical 
orbit of the satellite and compares a theoretical doopler 
curve obtained from this orbit with the actual dottier 
curve. The residuals are squared and summed and a new 
set of corrected variables is obtained. These allow the 
computation of a new theoretical orbit and a new theoretical 
doppler curve to be again compared with the observed data. 
This process is repeated until the sum of the scuares of 
the residuals is a minimum. At this point it is assumed 
that the best determination of the orbit has been made. 

In this way orbits of high precision can be determined 
and the future position of the satellite can be pr *dicte > 
with the required accuracy. It takes nearly two hours 
to update a single satellite orbit and prepare an injection 
message. With each satellite requiring this information 
approximately every twelve hours, the computer will be used 
to capacity. The message prepared by the compute!' contains 
^4,960 bits of oroital and time information which comprise 
the message proper and 1^0,000 bits of information which 
control the equipment at the injection facility during a 
satellite pass. 
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j,l though TRADAT .speed and accuracy are particularly 
-dapced to communications between the tracking stations 
and the operations/computer center, it still does rot 
provide the speed, and accuracy level reouired to handle 
the injection messages which are sent from the operations/ 
corrmuter carter to the injection facility. These must be 
100 per cent e^or free. High-speed, eouipment with a data 
rate more than twice that of TRADAT is used for the sole 
purpose of sending ii.jection messages, which, from the 
time they are fed into the equipment, are entirely automated 
until they are transmitted by the injection facility 
antenna to the satellite. 

4.5 THU INJECTION FACILITY. - - The operational 
system employs two injection facilities, one in California 
and one in Minnesota. Each injection facility contains 
a 60 foot oar fool ic directional antenna in conjunction 
with a 10 kilowatt transmitter . A small computer similar 
to that used in the AI\ T /BRN-3 submarine navigation eouipment 
controls the operation. The injection raeir ge from the 
operations/computer cento--'* is received at c,e injection 
facility and stored, Then the satellite is in range of the 
injection facility the message proper containing bits of 
information in combinations of binary ones and zeros is 
converted to waveform to modulate the carrier wave of the 
transmitter. The message is transmitted into the satellite 
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in 15 seconds. Since tho satellite and the earth are in 
motion during transmission, the distance from, the transmitter 
to the satellite varies constantly. Therefore, the time 
reouired for the signal to reach the satellite must be 
con-routed in advance and the bit rate must be adjusted to 
ensure that the satellite receives them at a even rate. 

The bit rate is changed 32 times during the 15 second 
transmission. The injection facility also reads the message 
as retransmitted from the satellite to determine che success 
of the injection attempt. If necessary, the injection 
facility recycles to do the whole operaoicn over again, 
complete cycle for each injection requires two minuses. 

The equipments of the injection facility consist of 
commercially available components wherever possible. Ir. 

fact the transmitter is a modified commercial television 

[&] 

transmitter . 

4.6 NAVIGATION TO, TT IPi/WT T T . - - Two types of shipboard 

navigation equipments have been developed, one for submarines 
and the other for surface ships. Tne AN / TV - - 3 was develcped 
by West inghouse and is designed specifically for submarine 
use. Ip incorporates many features adapted to submarine 
environm-nt, and certain pieces of interface equipment to 
permit operation in conjunction with other navigational 
equipments presently used aboard Fleet Ballistic r issile 
submarines. The AN/BRN-3 is completely automatic in 
operation. It computes it's own alerts of wr.vn unvigationul 



sat. ell it es will bo within radio range of the submarine 
position. It then computes the satellite position for 
specific intervals during a pass based upon knowledge it 
already has of the satellite orbit. Using an assumed 
submarine position it precomputes the expected doppier 
curve. When the actual doppier curve is received the 
AN/BRK-5 adjusts and readjusts the assumed position until 
the closest possible matching of the two curves is achieved 
Since the assumed position is supplied by an Inertial or 
some other high quality positioning system, the curve- 
fitting process takes place almost in real time. The 
position at which the two curve? fit most closely together 
is taken to be the submarine's Position. This fix is 
transmitted to the other navigation equipment on board, as 
necessary. In addition, the eouipment performs its own 
readiness checks, diagnoses any trouble, and prints out in 
english what part to replace. 

The eouipment developed for surface ships is designate 
the AN/ SIN- 9 . It is far less complex than the submarine 
equipment since it does not have the environmental and 
interface problems associated with undersea operations. 

The AN /3 RN - 9 receives the signal from tne satellite and 
automatically computes tne ship position. It uses the 
orbital parameters provided by the satellite to compute 
the part of the orbit corresponding to the doopler curve 
received. From this it computes the satellite coordinates 



65 



at the observation time points. 



plus the 
2 minute 

dY , and 

n 7 

position 

latitude 



Y.'ith this inf ormation, 

integrated doppler cycle count for at least three 

time intervals, it commutes the corrections dX , 

7 n 7 

d Z to the assumed position coordinates. The new 
n 

coordinates are then used to calculate the navigator 
and longitude. 



66 



5. SYSTEM 



ACCURACY 



5.1 GENERAL . - - In the course of the development 

of the Navy Navigation Satellite system, individual error 
components were analyzed in order to determine and improve 
those components which most affect the system. Theoretical 
values based on assumed idealized conditions were determined 
for most components. These values can be expected to vary 
under operating conditions. Specific values pertaining 
to the accuracy of fix determinations based on the use of 
the prototype and operational systems are not available 
due to security restrictions. However, several fleet 
units have made extensive use of the operational system 
and fix accuracies of i 0.1 mile or better have been made 
public . 

From the navigator's point of view the primary error 
components of the operational system are: 

1. Ionospheric refraction error. 

2. Frequency error. 

3. Time error. 

4. Error in the determination of the velocity 
of electromagnetic propagation. 

5. Navigator instrument error. 

6. Navigator velocity er^or. 

7. Satellite coordinate error. 

A short discussion of these error components is given in 
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the following sections. It should bo noted that duo to 
the nature of the system, the satellite coordinate error 
can be broken down into individual error components similar 
to those encountered in fix determination. The principal 
difference is in the high order gravity coefficients 
required for orbit determination and prediction. 

5.2 IONOSPHERIC REFRACTION. - - In Section 3.2 it 

was shown that a first order correction for ionospheric 

refraction can be determined by the transmission of two 

coherent frequencies from the satellite. The resulting 

doppler shift still contains higher order effects of the 

series (l/f ). In order to determine whether contributions 
s 

higher than the first order are significant, the theoretical 
effect of refraction on the doppler shift was considered 
in Reference 9. Table 5.1 summarizes the results pertaining 
to upper bounds on the various contributions to the refracted 
doppler shift for various frequencies within the transmission 

DO 

range. Except for the vacuum term, typical values are 
probably a factor of three to ten smaller. The first order 
term is proportional to the time rate of change of the 
total number of electrons in a unit cross section tube 
extending from the observation point to the satellite. This 
is the term that is eliminated when two frequency doppler 
data is used. The second order term considers Faraday 
rotation and its sign denends on whether the antenna is left 
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or right circulary polarizod. The third order correction 
is composed of four terms. The first term arises from 
the fourth order term in the expansion of the refractive 
index series. The second term represents the difference 
between the optical path and the geometric path. The Inst 
two terms represent modifications to the solution o^ the 
wave equation using the geometrical optics approximation. 
These terms depend upon various derivatives of the electron 
density. Unless the satellite frequencies are below VHP 
the two terms not obtainable from the geometrical optics 
approximation are negligible. 

Table 5.2 presents the estimated maxima of the second 
and third order contributions to the vacuum doppler shift 
when dual frecuency data are combined to eliminate the 

\3 I 

first order refraction term. Three typical freauency 
combinations ace given and the table entries have been 
converted from frecuency shift to errors in the range rate 
in meters per second. 

Examination of the tables indicates that if the lower 
of the two frequencies is not less than 100 mc/s and the 
ionosphere is not extremely disturbed, the second and third 
order refraction contribut ions should be negligible. 

5.3 FREQUENCY AND TIKE ERRORS . - - The navigator 

must have an accurate clock and a reference frecuency in 
order to determine his position from doppler observations. 
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In addition he must know the freouency transmitted by the 
satellite and the satellite position as a function time. 

In Chapter 3 it was shown that to the required accur cy the 
frequency of the reference oscillator is assumed unknown 
and the position determination computations were arranged 
so as to eliminate the need for this value. The navigator 
is able to recover time from each satellite observation to 
an accuracy of better than one millisecond. This time is 
maintained by a clock derived from the reference oscillator 
which has a drift rate of less than 1 part in 10^ per day. 
Therefore, the time and freouency errors in the Navy 
Navigation Satellite System both derive from the drift of the 
satellite oscillator. A drifting oscillator gives an along 
track navigation error because the frequency is in error 
and causes a time error in the clock that' it is driving. 

This time error ."Iso gives an along, track error. Navigation 
errors of about 1 kilometer result from an error of 1 part 
in lCr in the satellite frequency. However, satellite 
freouency can be determined with an accuracy of at least 

1 part in 10^ and predicted with an accuracy of at least 

2 parts in 10^ per day. Since the satellite oscillator 
can be tuned by command from an injection station twice 
per day, it’s freouency can bo keDt to a predetermined 
value within an accuracy of a least 2 parts in 10^. This 

D°] 

causes a navigation error of approximately 20 meters. 
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6.4 VELOCIT v OF ELECTROI^GNET IC PROPAGATION. 



The 



value for the velocity of electromagnetic energy propagation 

is taken to he a constant in calculation® for nosition 

determination using the Navy Navigation Satellite System. 

The value is probably known to an accuracy of 1 part in 

3,000,000. This world oermit position determination to an 

[hi 

accuracy of about two meters. ' Future improvements in this 
constant can be accomodated by modification of the system 
based on the current accepted value. 

5.5 NAVIGATION INSTRUMENT ERupR. - - In order to test 
the contribution that the navigation equipment makes to the 
navigation error, observations were made on experimental 
satellites transmitting on four frequencies. Two navigation 
instruments at the same site made simultaneous observations 
on the satellite. One in^trTinent used two of the satellite 
freauenci es and the other used the remaining two frequencies. 
The difference between the navigated oo®itions is then the 
result oT two independent measurements of position. The 
result of one set of tests obtained with the Transit 2A 
satellite gave a mean difference in position of 50 meters 

QQ 

in one coordinate and 30 meters in the other. It is 
expected that instrumentation has improved since this time. 

5.6 NAVIGATOR VELOCITY ERROR. - - Since the navigator 
must u®e a measured value of his velocity to determine his 
position at the times of the satellite observations , an 
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error in this measurement will cause an er^or in the fix. A 
ship acted on by unknown current, windage drift, and sea 
condition effects amounting to a one knot velocity error ^.n 
the northerly direction has an estimated position error of 



causes a negligible positional error at mid-latitudes since 
the satellite is in a polar orbit. It was determined to at 
least three place accuracy that the position errors varied 

GO 

linearly with the velocity errors. This does not appear 
to be a severe limitation for surface ships and submarines 
where accurate pit-log or inertial system velocities are 
normally available. 



preceeding error components which are normally negligible, 
the only remaining source of system er”or is in the satellit 
coordinates transmitted to the user. Any error in these 
coordinates is reflected almost exactly one for one in the 
accuracy of the navigation fix. Tests on the size of the 
errors in the predicted orbits of the Navy Navigation 
Satellite System were conducted by computing and predicting 
the orbits at four tracking stations. These predictions 
were made using geodetic parameters in use at the Naval 
Weapons Laboratory during four periods of time through 1964. 
The predicted orbit was compared with an orbit based upon 
independent doppler observations made by world-wide tracking 
stations during the prediction interval and using the latest 



0.3 nautical miles 




A one knot easterly velocity error 



5.7 SAT -ILL IT'S COORWATR RR^OR. 



Apart from the 
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available gravity coefficients. The resulting twelve 
hour prediction errors for random sets of observations 
made on polar orbiting satellites show a decrease in 
size from over 1000 meters in 1960 to the value of 50 



meters in 1965. 



DO 
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5.8 PROPOSE 4RU0R ANALYSIS . - - A high speed comouter 
may be easily programed to analyze the effect of individual 
and/or collective error components on a given satellite 
position determination. The basic inputs of this program 
would consist of the following parameters: 




2. iAx, iAy, iZz 

5. Xf, Y l, Z l 

4. f r 

5 - f s 

6 . c 

7. AT 

The satellite coordinates as a function of time in an 
earth centered rectangular coordinate system are available 
from a Geodetic Science Computer Library Program at the 
Ohio State University or from some other similar program. 
Using satellite coordinates bracketing three or more time 
intervals and choosing the above parameters which are assumed 
errorless, we cun determine the expected integrated cycle 
count from (3.7). Assuming errorless instrumental and 
environmental conditions the calculated cycle counts and 
the above oarameters may be considered as an ideal errorless 
position determination based on satellite observations. 

This can be verified by substitution of these values into 
the equations developed in Chapter 3. The position 
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coordinates produced by calculation will be equal to the 

initial navigator coordinates since the corrections dX n , 

dY . and dYf will be zero in our assumed errorless 
n n 

condition. We have thus created an ideal errorless 
position determination which can be used as a standard 
in an error analysis. By random and/or systematic 
alteration of the assumed errorless input parameters we 
can determine their individual and/or combined effect on 
the final position determination. A study could also be 
made to determine the effect of decreasing the time interval 
and using more observations in a least squares solution. 
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summary 



6 . 



The adoption of space technology provides a new aporoach 
to the earth-bound problem of marine navigation. The Navy 
Navigation Satellite System is based on the doppler shift 
in a radio frequency transmitted from a near-earth satellite. 
By exact measurement of the doppler shift it is possible to 
determine the location of the receiver on the surface of the 
assumed figure of the earth. The accuracy obtained by using 
this technique is possible because the satellite is very 
precisely restricted to an orbit defined by physical laws 
of motion and the quantities observed, time and freouency, 
can be measured to an accuracy of one part in 10^^. Of all 
the possible satellite paths permitted there is only one 
which results in a particular curve of doppler shift. 

Two basic problems are inherent in this system. One 
is that the ionosphere bends radio waves and thus gives a 
false position of the satellite. This has been overcome 
by the use of two satellite transmission fr>eouencies which 
allow a first order correction to be computed. The second 
basic problem is the determination of the gravity field of 
the earth. This affects the determination and prediction 
of satellite orbits. This problem has been reduced to an 
acceptable level by using satellite orbital data to calculate 
new gravity coefficients. 

The use of a navigational satellite system based on 
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the dopoler shift of VHP radio waves appeals to nrovide 
the beat present day solution to the military navigation 
requirements. World -wide coverage is obtained by use of 
polar orbits. The system is inherently all weather and 
immune to interference since line-of-sight radio frequencies 
are used for both tracking and navigation. lince the 
satellites transmit their ephemerides the navigator need 
neither interrogate the satellite nor receive orbital and 
time information by other communication links. The system 
is also reliable because a temporary interruption of 
reception during a pass does not preclude the use of 
data that are received to obtain an accurate navigational 
fix. It is assumed that the system meets the military 
accuracy reauir ements . 
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